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Abstract. In the past years we have made great efforts to reduce the sta- 
tistical and systematic uncertainties in stellar parameter and chemical abun- 
dance determinations of early B-type stars. Both the construction of robust 
model atoms for non-LTE line-formation calculations and a novel self-consistent 
spectral analysis methodology were decisive to achieve results of unprecedented 
precision. They were extensively tested and applied to high-quality spectra of 
stars from OB associations and the field in the solar neighbourhood, covering 
a broad parameter range. Initially, most lines of hydrogen, helium and carbon 
in the optical/near-IR spectral range were reproduced simultaneously in a con- 
sistent way for the first time, improving drastically on the accuracy of results 
in published work. By taking additional ionization equilibria of oxygen, neon, 
silicon and iron into account, uncertainties as low as ~1% in effective tempera- 
ture, ~10% in surface gravity and ~20% in elemental abundances are achieved 
- compared to ~5-10%, ~25% and a factor ~2-3 using standard methods. 

Several sources of systematic errors have been identified when comparing 
our methods for early B-type stars with standard techniques used in the nineties 
and also recently (e.g. VLT-FLAMES survey of massive stars). Improvements 
in automatic analyses are strongly recommended for meaningful comparisons of 
spectroscopic stellar parameters and chemical abundances ('observational con- 
strains') with predictions of stellar and galactochemical evolution models. 



1. Introduction 

Normal unevolved early OB-type stars of ~8-20 M are the objects with the 
simplest photospheric physics among the massive stars. They are unaffected 
by e.g. strong stellar winds like the hotter and more luminous stars or by 
convection and chromospheres like the cool supergiants. However, their spectral 
analysis turned out to provide inconclusive results in the past decades, i.e. too 
large uncertainties in basic stellar parameters and an overall enormous range in 
derived elemental abundances, posing a challenge to predictions of stellar and 
Galactochemical evolution models (see review by Przybilla 2008). 

In order to improve the quantitative analysis of these stars we have exhaus- 
tively updated the spectral modelling by constructructing robust model atoms 
for non-LTE line-formation calculations. In parallel, we have implemented a 
powerful self-consistent analysis technique, which brings numerous spectroscopic 
parameter and abundance indicators into agreement simultaneously. 

Our efforts have provided highly-promising results so far, i.e. a drastic 
reduction of statistical and systematic uncertainties in stellar parameters and 
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chemical abundances (Nieva & Przybilla 2007, 2008; NP07/08). As a first ap- 
plication, stars from OB associations and the field in the solar neighbourhood 
covering a broad parameter range were analysed. The sample turned out to be 
chemically homogeneous on the ~10% level (Przybilla, Nieva & Butler 2008, 
PNB08), corroborating earlier findings from analyses of the ISM gas-phase. The 
data are also consistent with published Orion nebula abundances. The results 
have an immediate impact on several fields of contemporary astrophysics like 
stellar (see Przybilla, Firnstein & Nieva, these proceedings) and Galactic chem- 
ical evolution models and the dust-phase composition of the local ISM. They 
provide an independent view on the discussion of photospheric solar abundances 
and helioseismic constraints on the solar interior model, and they define the ini- 
tial chemical composition for models of star and planet formation in the solar 
neighbourhood. In addition, a few B-type hyper-velocity stars were analysed 
using this technique, providing valuable constraints on their nature and their 
ejection mechanisms (Przybilla et al. 2008a, b). 

Despite this kind of star has relatively simple photospheres compared to 
other objects, the spectral analysis is still sensitive to many potential systematic 
effects which are usually underestimated. We discuss the most common sources 
of systematic error that have to be avoided when high precision/accuracy in spec- 
tral analyses is desired. Observational constraints as obtained from automatic 
spectral analyses of large star samples like the VLT-FLAMES survey of massive 
stars may benefit significantly from a proper elimination of these systematics. 

2. Model calculations and new spectral analysis methodology 

A hybrid approach is used for the non-LTE line- formation computations. These 
are based on line-blanketed plane-parallel, homogeneous and hydrostatic LTE 
model atmospheres calculated with Atlas9. Non-LTE synthetic spectra are 
computed with recent versions of Detail and Surface. These codes solve the 
coupled radiative transfer and statistical equilibrium equations and compute 
synthetic spectra using refined line-broadening data, respectively. The hybrid 
non-LTE approach is consistent with full non-LTE calculations (NP07) but faster 
and it also allows comprehensive model atoms based on critically selected atomic 
data to be employed in the non-LTE line- formation computations. 

The new spectral analysis was originally based on a self-consistent and si- 
multaneous reproduction of almost all hydrogen, helium and carbon lines in 
the optical/near-IR spectra, matching multiple ionization equilibria (Hel/ll, 
Cii/iii/iv), see Nieva k Przybilla (2006), NP07 and NP08 for details. The 
method was further extended by consideration of additional ionization equilib- 
ria (Oi/ii, Nei/n, Si ii/iii/iv, Feii/m, PNB08). This allows unprecedently 
accurate stellar parameters and elemental abundances to be derived, with un- 
certainties as low as ~1% in effective temperature T e s , ~10% in surface gravity 
g and ~20% in elemental abundances. Significant improvements on results from 
previous studies are thus achieved, which typically give uncertainties of ~5-10%, 
~25% and a factor ~2-3 for these quantities. Moreover, the successful imple- 
mentation of the new method required an identification of sources of systematic 
uncertainties in standard spectral analyses and allowed in most cases for their 
quantification, which we discuss in the following. 
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3. Reducing systematic uncertainties 

Every step in a quantitative spectral analysis is susceptible to systematic un- 
certainties. Here we briefly list the most common sources of systematics that 
affect the final error in stellar parameters and elemental abundances of normal 
unevolved early B-type stars. When possible, recipies are given how to prevent 
them. Additional systematics may arise from further complications like mag- 
netic fields, but this is beyond the present scope. 

Single or double? Not only spectroscopic but also close visual binaries are 
expected to be observed in dense fields. Standard analyses for single stars ap- 
plied to a spectrum contaminated with light of a second star can give erroneous 
results throughout. Inspection of H or He lines can help to identify asymmetries 
due to a companion before earring out an automatised quantitative analysis. 

Quality of spectra. Continuum normalization, local continuum definition and 
low S/N are important sources of systematics. E.g., spectra of S/N~50 challenge 
the definition of the local continuum of spectral lines. The abundance determi- 
nation in stars rotating at intermediate velocities (usini ~ 50kms _1 ) is limited 
to 0.2-0.3 dex in accuracy at this S/N. Fast-rotating stars (vs'mi > 150 km s -1 ) 
or low-resolution spectra impose even more complications because metal line 
blends lower the real continuum, hence the abundances can be systematically 
underestimated and the accuracy is limited to 0.3-0.4 dex (see Korn et al. 2005). 

Model atmospheres and line formation. Atmospheric structures computed 
with full non-LTE or hybrid non-LTE methods are equivalent for unevolved B 
stars (NP07), when the abundances used for opacity calculations are the same - 
note that 'solar' abundances have changed over time. Line-blanketing effects on 
the atmosphere impact the line- formation calculations, introducing dependen- 
cies on metallicity and microturbulence. In a similar way, line blocking affects 
the strength of synthetic lines by modifying radiative rates. Both effects need 
to be accounted for in a consistent way. 

Model atoms. Most model atoms for non-LTE calculations are based on input 
atomic data (ab-initio and approximation data) as available in the early nineties. 
It is worthwhile to check for improvements on the modelling whenever new data 
becomes available. E.g., for C different model atoms yield discrepancies in abun- 
dances up to 0.8 dex for some lines while no discrepancies are found for others, 
see Sect. [3~Ll and NP08. Other elements also show a similar behaviour. 

Effective temperatures. T c fj estimated from photometry can differ from those 
determined with the self-consistent spectroscopic method by more than 10% 
(NP08). Spectroscopic determinations via ionization equilibria are a powerful 
technique only when the model atoms are reliable. In addition, one ionization 
equilibrium alone does not provide accurate constraints because of dependencies 
on other variables like microturbulence (see below, Fig. [2] and consequences for 
C abundances in Table [1]). Therefore, only the simultaneous use of multiple 
ionization equilibria provides a reliable and parameter-free T e g- determination. 

Surface gravity. The common use of only one Balmer line as log g indicator 
does not allow for consistency checks. Instead, most Balmer lines and metal ion- 
ization equilibria should be considered for an accurate determination. Moreover, 
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Figure 1. Left: sensitivity of two Cn lines to use of photoionization cross- 
sections from different ab-initio calculations in a model spectrum for a B2 V 
star. Right: differences in C abundance as derived from the two lines in 6 
stars as a function of T c s . From NP08, see text for details. 



neglecting non-LTE effects also leads to systematic errors: increasing with T e g 
up to ~0.2dex in log g around 35 000 K, see NP07/08 and Sect. 13.1.1 Ionization 
equilibria can also be a powerful tool for the log g determination. 

Microturbulence. This quantity is generally derived from selected lines of 
one oxygen or one silicon ion. However, the microturbulent velocity £ should 
adopt the same value for all metals. In particular, cross-checks including dif- 
ferent species are mandatory when only a few lines are measurable (like in fast 
rotators) in order to avoid large uncertainties. Microturbulent velocities often 
come close to or exceed the sound speed in many previous studies. In contrast, 
physically more reasonable lower values are found in our work, increasing from 
~2-4kms _1 in dwarfs to ~5-8kms _1 in giants. Consequences of overestimated 
£ for the T e g and abundance determination are discussed in Sect. 13.1.1 

Spectral line selection. Abundances may depend on the selection of lines 
used in the analysis when the model atoms are not comprehensive: some multi- 
plets may indicate systematically different abundances than others. Moreover, 
some lines are observable only in cooler or in hotter stars. Which lines should 
one choose for the analysis? All possible observed lines for each star should be 
taken into account in the optimum case (see NP08 and PNB08). 

Non-LTE corrections. Non-LTE effects cannot be easily predicted. They 
affect different lines of the same star in different ways. Non-LTE line strength- 
ening or weakening can occur, and non-LTE effects are not restricted to stronger 
lines alone. For different plasma conditions the non-LTE effects change. Hence, 
adding or subtracting fixed 'non-LTE abundance corrections' to LTE results may 
increase the systematics. 

Macroturbulence. Macroturbulence is not considered in many studies. This 
is important for proper v sini determinations in apparently slow-rotating objects. 
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Figure 2. Effective temperature determination via Si ionization equilibrium 
for the LMC star NGC2004-D15. The Sim lines are strong, hence the de- 
rived Si in abundance depends on the microturbulence £. Different T e g is 
derived from Si II / III and Sim/iv for ill-chosen £, e.g. an overestimation of 
£ by 7kms _1 gives AT e g «2000K. The Sill/m/iv ionization equilibrium is 
established only for the correct £ (from Nieva 2007). 



3.1. Systematics from atomic data and atmospheric parameters 

Here we provide a few quantitative examples of systematic effects on abundance 
analyses to be expected due to use of different input atomic data for non-LTE 
calculations and due to atmospheric parameter variations. 

Fig. Q] shows discrepant results for C n lines predicted by two different non- 
LTE model atoms. On the left panel the sensitivity of two lines to different 
ab-initio photoionization cross-sections is shown. Inaccurate atomic data will 
indicate largely underestimated abundances for the stronger line. This is quan- 
tified for stars within a broad parameter range (right panel). The discrepancies 
in abundances depend on the reliability of the model atom and on the plasma 
conditions of the star and amount up to 0.8 dex. In parallel, collisional excitation 
and ionization cross-sections and log gf values need to be reliable (NP08). 

Fig. [2] shows the dependency of T e g to the adopted value of microturbulence 
when only one ionization equilibrium of Si is adopted, i.e. Sin/in or Silll/iv. 
This will also affect the derivation of log g from the Balmer lines. This problem 
can be solved via use of multiple ionization equilibria, e.g. Siii/in/iv and when 
possible also considering other elements (as explained above). The goal should 
be deriving the same value of T e ff , log g and £ from all H, He and all metals. 

Systematic errors in C abundances from individual lines in a Bl III star due 
to variations in T e g, log g and £ have been quantified in Table [TJ The offsets in 
parameters are averaged discrepancies for studies using standard analysis tech- 
niques. Incorrect parameters prevent consistent Cn/in/iv ionization equilibria 
to be achieved, even when the model atom is highly reliable in the defined pa- 
rameter range. Parameters derived from C ionization equilibrium (NP08) are 
confirmed by other metals (PNB08) when reliable model atoms are used. 



4. Conclusions 

In the past years we have made great efforts to reduce uncertainties in quan- 
titative spectral analyses of early B-type stars. A self-consistent analysis, i.e. 
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Table 1. Systematic errors in C abundances of individual lines caused by 
atmospheric parameter variations and the assumption of LTE for the lino- 
formation calculations in the BIHI star HR3055 (T off = 31 200 ± 300 K; 
logg = 3.95 ±0.05; £ = 8±2km S -\ see NP08, PNB08) 



Ion 


A 


W x 


e(C)nlte 


AT eff 


A log g 


AC 


LTE 




(A) 


(mA) 




-2000 K 


+0.2 dex 


+5kms 1 




Cii 


4267.2 


113 


8.46 


-0.33 


-0.11 


-0.16 


-0.40 




5133.3 


19 


8.34 


-0.30 


-0.10 


0.00 


0.00 




5143.4 


10 


8.44 


-0.40 


-0.05 


0.00 


0.00 




5145.2 


19 


8.36 


-0.32 


-0.09 


-0.02 


0.00 




5151.1 


11 


8.34 


-0.30 


-0.08 


0.00 


0.00 




5662.5 


4 


8.37 


-0.33 


-0.13 


0.00 


0.00 




6578.0 


66 


8.24 


-0.40 


-0.15 


-0.10 


0.00 




6582.9 


37 


8.29 


-0.30 


+0.02 


+0.05 


+0.05 


Cm 


4056.1 


45 


8.28 


+0.21 


+0.06 


-0.04 


+0.08 




4162.9 


58 


8.29 


+0.28 


+0.09 


-0.03 


+0.25 




4186.9 


92 


8.34 


+0.35 


+0.15 


-0.08 


+0.07 




4663.5 


18 


8.27 


+0.22 


+0.07 


-0.03 


+0.22 




4665.9 


50 


8.24 


+0.26 


+0.08 


-0.08 


+0.35 




5272.5 


14 


8.28 


+0.16 


+0.01 


0.00 


0.00 


Civ 


5801.3 


53 


8.45 


+1.06 


+0.46 


-0.03 


+0.39 




5811.9 


34 


8.45 


+1.06 


+0.46 


-0.03 


+0.39 



account of all spectroscopic indicators - Balmer and helium lines and multi- 
ple metal ionization equilibria - throughout the optical and near-IR, resulted 
in drastically reduced systematic effects in the atmospheric parameter and ele- 
mental abundance determination. A large number of potetial systematic errors 
was identified when comparing our new models and self-consistent method with 
standard techniques. One should keep in mind that statistics does by no mean 
reduce systematic errors. We conclude that careful improvements in automatic 
spectral analysis routines should be implemented before being applied to large 
samples of stellar spectra. This will prevent unnecessary systematic bias in 
stellar parameters and abundance determinations, inconclusive results and mis- 
interpretations of the studies by theoreticians. Comparisons of 'observational 
constraints' and theoretical predictions are only meaningful when the former 
are unbiased by systematic error. 
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